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ABSTRACT 

We study and constraint Mass- Varying Neutrino models using present and future available data. In 
these models, dark energy is a self-interacting scalar field directly coupled to neutrinos. We investigate 
two different potentials and both positive and negative coupling parameter /3. This corresponds to 
increasing or decreasing neutrino mass, respectively. We explore couplings up to < 5. In the 
case of the exponential potential, we find upper limits on < 0.004 at 2-a level. In the case of the 
inverse power law potential the null coupling can be excluded with more than 2-a significance, the 
limits on the coupling being /3 > 3 for the increasing neutrino mass and /3 < —1.5 for the decreasing 
mass case. This is a clear sign of a preference for higer couplings. When including a prior on the 
neutrino mass today the upper limits on the coupling become < 3 at 2-(j level for the exponential 
potential. Finally, we present Fisher forecast using the tomographic weak lensing from the Euclid- 
like experiment, also in combination with the CMB temperature and polarization spectra from the 
Planck-like mission. If considered alone, lensing data is very efficient in constraining Wi/, giving a 
signal of a non-null neutrino mass with high significance. There is, however, a strong degeneracy in 
the P-uJu plane. When the two data sets are combined, the latter degeneracy remains, but the zero /3 
value can be excluded at more than 2-a. 



I. INTRODUCTION 

Present days astronomical observations strongly indi- 
cate there are two unknown components in the universe 
matter-energy budg et. These are dark matter (DM) and 
dark energy (DEldRiess et all 119981 : IPerlmutter et all 
[1991 [Kowilski et al.ll20'08i r The cosmological constant 
A is one of the simplest candidates for the nature of dark 
energy. Although so far consistent with al l major ob- 
servationa l probes (|Wood-Vasev et al.ll2007t IReid et al.l 
[2010: Komatsu et al.ll2011l : ISuUivan et al.ll201lD . it faces 
two major puzzles: 1) Why such a tiny value (the cos- 
mological constant problem)? 2) Why A has become 
important only recently (the coincidence problem) . Sev- 
eral possible solutions to these issues have been put for- 
ward. In p articular, by identifying dark energy w ith a 
scalar field (jWetterichI 11981 [Peebles fc Ratral[l988l ) and 
to cou ple it to either dark matter, or baryons or neu- 
trinos (lWetterichlll995l:IAmendolall2000l:IAmendola et all 
120111: [La Vacca et al.ll2009t iKristiansen et al.ll2oTofr 

In this work, we study models where the dark 
energy field is coupled to neutrinos. These mod- 
els are called Mass Varying Neutrino (MaVaN) 
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indeed be solved for a mass growing neutrino with a cou- 
pling to quintessence somewhat larger than gravity. 

Several authors have shown that important cos- 
mological effects ca n ap p ear w ithin this class 
of models (Biaeld e et al.[ [2008[: [ feookficld et al 
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I2011al: [BaldH [20lll : lAvaita et all [20ll . Specially, 
when the neutrino mass is sufficiently big for neu- 
trinos to be non-relativistic. When this happens, 
neutrinos feel the presence of the fifth force and can 
collapse i nto nonlinear struc t ures, w hich are stable an d 
bounded [Winte rgcrst et al.' (2010): iWettericlj ([20071); 
[Bernardini fcBertolami (2009, 2008 ). It has been shown 
that these neutrino lumps form at redshift Zni ~ 1 — 2, 
when the neut r ino fluctua t ions become nonlinear 
([Wetterich[ [2007t [Mota et all [2008t [Wintergerst et all 
I2O1O0 . These effects are mainly present in the strong 
coupling regim e where a simple line ar approach may 
be insufficient ([Pettorino et al.l [20101 ). Here we limit 
our analysis to the low coupling regime in which such 
instabilities are absent. 

In this work we further investigate the cosmological sig- 
natures of MaVaN models and put constraints on them 
using the most updated observational data from the Cos- 
mic Microwave Background (CMB) Radiation tempera- 
ture anisotropics spectra, from the Large Scale Struc- 
ture (LSS), and the Supernovae Type la (SNIa) lumi- 
nosity distance. We start by investigating two differ- 
ent scalar field potentials and both positive and nega- 
tive coupling parameter. This corresponds to increas- 
ing or decreasing neutrino mass, respectively. Finally, in 
the last section we present Fisher forecast using the to- 
mographic weak lens ing from an Euclid-like experiment 
([Laureiis et al.|[2"011f ). also in combination with the CMB 
tempera ture and polarization spectra fro m a Planck-like 
mission ([The Planck Collaboratioii|[2006[ ). 

2. THEORY 

2.1. The Cosmological Background Evolution 

MaVaN models involve a coupling between the DE 
scalar field and massive neutrinos. In a flat, homoge- 
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neous, Friedmann-Robertson- Walker universe with line- 
element 



(1) 



the Friedmann equation describes the universe expan- 
sion: 



(2) 



In this equation, Pa{T) is the energy density of the in- 
dividual components a, including CDM, DE, neutrinos, 
baryons and radiation. The prime refers to the derivative 
with respect to conformal time t. 

The coupling of DE to the neutrinos results in the neu- 
trino mass becoming a function of the scalar field, whose 
evolution is described by the coupling expression 
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We choose to take a constant /3 such that: 



(3) 



(4) 



where is a constant. An additional attractive force 
between neutrinos of strength 2/3^ is mediated by the 
scalar field exchange. For /3 ^ 1 this corresponds to a 
strength comparable to gravity. 

The existence of the coupling between DE and massive 
neutrinos can be formalized, assuming that the sum of 
the energy momentum tensors for the two species is con- 
served, but not the separate parts. We neglect a possible 
scalar field coupling to CDM, so that p'^ = —3Tipc- This 
hypothesis requires DE and mass varying neutrinos obey 
the coupled conservation equations: 

p'^ = -3nil + w^)p^+l3(j)'{l-3w,)p, (5) 

pI ^ -3H{1 + 'Wiy)p^ - [3(I)'{1 - 3wv)pi, (6) 

where, the energy density stored in the neutrinos is given 

by 

P. = \ I q'dqdnefoiq), (7) 



and the pressure by 



1 

3^ 



q^dq dnfo{q)- 



(8) 



where fo{q) is the usual unperturbed background neu- 
trino Fermi-Dirac distribution function 



/o(e) 



1 



(9) 



and = q^ +TO^((/))a^ {q denotes the comoving momen- 
tum). As usual, hp and ks stand for the number 
of spin degrees of freedom, Planck's constant and Boltz- 
mann's constant respectively. In the following we will 
assume that the neutrinos decouple whilst they are still 
relativistic, and therefore the phase-space density only 
depends upon the comoving momentum. 

The energy density and pressure of the quintessence 
field are given by 

= + = g - (10) 
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Fig. 1. — Energy density evolution of di ffere nt components, as 
specified in the legend, for the potential I I12I I. with 13 = —1.9, 
a = 0.021, uj^ = 0.0010, ujcdm = 0.1097, ujt = 0.0224, ho = 70. 

Taking into account the energy conservation of the cou- 
pled neutrino-dark energy system, one can immediately 
find that the evolution of the scalar field is described by 
a modified Klein-Gordon equation: 



2-H0' 



,dV 
d0 



a'pip, - 3p,) . (11) 



This equation contains an extra source term with respect 
to the uncoupled case, which accounts for the energy 
exchange between the neutrinos and the scalar field. 

Here we investigate two different expressions for the 
DE potential: an expo nential potential (fWetterich. 19881 : 
iFerreira fc Jwci[l998l) : 



and an inve rse 
iBinetruvl [19991 ) 



y(0) oc e""'^ (12) 
power- law (|Steinhardt et al.l 119991: 



V(0) cx (A//0)2 



(13) 

Both a in (fT2| and M in (|13p are constant. In particular 
a is a free parameter which determines the slope of the 
potential and thus the DE fraction at early times. On 
the other hand, M defines the potential energy scale and, 
together with in Eq. Ul is used to set the convergence 
of the model to the actual cosmology. 

Once the potential (fT2|) or ([T3|) is given, the evolution 
equations can be numerically solved. In the very early 
universe, neutrinos are still relativistic and almost mass- 
less, with = Pi^/S and the coupling term in eqs.([51), 
(HI), ini) vanishes. 

However, as soon as neutrinos become nonrelativistic, 
the coupling term ~ becomes significantly different 
from zero, affecting also the evolution of the field (f>. As 
shown in Fig. [U pi, and p^ change behavior for z < 6: 
the value of the scalar field stays 'almost' constant and 
the frozen scalar field potential mimics a cosmological 
constant. According to this model, at the present time 
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neutrinos are still subdominant with respect to CDM, 
though in the future they will take the lead. For our 
choice of the parameters neutrino pressure terms may be 
safely neglected for redshifts Znr < 4; before that redshift 
neutrinos free stream as usual relativistic particles. 

2.2. Perturbed Equations 

In order to compute the evolution of cosmological per- 
turbations in our model, we work in the synchronous 
gauge, taking the line element to be 



ds^ 



hij) dx^dx-' . 



(14) 

the perturbed Klein-G ordon equation is given by 
()Brookfield et al.ll2006bD : 



2HS(t> 
dlnm 



(15) 



d(t> 



-{6p„ - 3dpu) + 
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For the neutrinos we use the perturbed part of the 
energy momentum conservation equation for the coupled 
neutrinos 

_ d\nm 

to calculate the evolution equations for the neutrino per- 
turbations (T*^ stands for the trace of the neutrino en- 
ergy momentum tensor). Taking 7 = we derive the 
equation governing the evolution of the neutrino den- 
sity contrast , by = whilst taking 7 = i (spatial in- 
dex) yields the velocity perturbation equation dy = ikiV^, 
with the coordinate velocity vl = dx"^ /dr: 



+ /3 (1 - 3wy) S<j) + '^^5(j> (1 - iwy) 



(17) 



1 + 



1 + Wy 



+ 1^1 — '^-^eSct)- pil-iwy^e^-k^cfy. (18) 

The variable represents the neutrino anisotropic 
stress. It is the presence of the additional coupling terms 
in these expressions for the growth of the neutrino den- 
sity and velocity perturbations, as well as the modifica- 
tions to the evolution of the cosmological background, 
which alters the behaviour of the neutrino perturbations 
in comparison with the standard uncoupled case. 

3. METHODS AND DATA 

Two are the methods we use to test MaVaN the- 
ory against cosmological data: the likelihood analysis 
through Markov-Chain Monte-Carlo (MCMC) technique 
and the Fisher information matrix. 

The CMB temperature and polarization anisotropy 
spectrum and matter powe r spectrum are ca lculated by 
suitably modifying CAMB ([Lewis et al.ll2000f) to contem- 
plate MaVaN's equations as described above. To ensure 



flat priors, unless other- 
choose to span only the 
< 5, since linear per- 



the accuracy of our calculations, we directly integrate 
the neutrino distribution function, rather than using the 
standard velocity weighted series approximation scheme. 

For the MCMC analysis, we use a mo dified version of 
the p ublicly available code COSMOMC (jLewis fc Bridld 
l2002f l to explore the parameter space. We consider the 
following basic set of parameters: 

{ Wfc , WCD A/ , ^^0 , Zre , ns , In As , , /?} , 

plus a when the exponential potential is involved. Here: 
uJb.cDM.v are the physical baryon, total dark matter and 
neutrino density parameters, uJb,CDM,v = ^b,CDM,vh'^ , 
where h is the dimensionless Hubble parameter Hq; z^e 
is the redshift of reionization; is the scalar spectral in- 
dex; As denotes the amplitude of the scalar fluctuations 
at a scale of k = 0.05 Mpc~^. The sum of v masses is di- 
rectly related to the neutrino density parameter through 
the relation A/^ = Ejti^ = • 93.5 eV, assuming 3 equal 
mass J/'s. 

All parameters are given 
wise stated explicitly. We 
low coupling regime up to 

turbations become unstable for high couplings, e.g. 
|/3| ~ 50. This is not a new issue in MaVaNs mod- 
els, where many other auth o rs have reported those 
instabilities (iMo ta et all l2008t iWintergerst et all 120101 : 
lAfshordi et al.l l2005: Fr anca et al.ll2009f l. However, there 
are several other claims in the literature that as the 
neutrino perturbations become n onlinear, those struc- 
ture virialise and become stable (jPettorino et all 120101 : 
iBaldi et al.ll2011bl ). Hence, the strongly coupled mass- 
varying neutrino models may not be ruled out by the 
large scale structure formation in the linear regime, such 
as the integrated Sachs- Wolfe effect or the matter power 
spectra. In spite of the many atter npts to understand and 
manage those "n on-linearities" (jPettorino et al.l 120101 : 
IBaldi et al.ll2011bD . a proper method of coping with these 
effects in a MCMC method to probe the parameter space 
of the models in the high /? regime is still lacking. There- 
fore, we leave the strong coupling regime for a future 
work. In our MCMC analysis we assume that the Uni- 
verse is spatially flat. 

With the aim of obtaining the best estimate of the 
cosmological parameters, we combine different CMB 
datasets with data from large scale structure and SNIa. 
In particular, for the CMB temperature anisotropy spec- 
tra we considered the 7 ye ar data from the WMA P satel- 
hte (WMAP7, see e.g. iKomatsu et~an (120111)) plu s 
higher m ultipole data from CB I (iSievers et al.l I2007T). 
ACBAR (jReichardt et al.ll2009f) . VSA (jPickinson et all 
l2004f l experiments. For the matter pow er spectrum we 
choose the results from the SDSS survey (jTegmark et al.l 
12006( 1. Constraints on the recent expansion history of 
the Universe are given b y the SNIa observation s from 
the Union2 compilation (jAmanullah et all I2010D . We 
also apply addit ional priors on t he Hubble parameter 
Ho = 74.2 ± 3.6 (|Riess et al.ll20M ). 

In addition to the likelihood analysis, we tested the 
power of future CMB and weak lensing data in con- 
strai ning MaVaN's parameters using t he Fisher formal- 
ism (|FisheHll935HTegmark et al.|[l99l . 

Fisher information measures the amount of informa- 
tion that observable random variables X — (Xi, . . . X„) 
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provide about unknown parameters ~ . . . 6^) upon 
which the probabiUty of X depends. Let C{'K.\6) denote 
the hkelihood function, yielding the probabihty that a 
model defined by the parameters 0, gives the set of ob- 
servables X. In our case X^ are CMB anisotropy and 
weak lensing spectra. Then the Fisher information ma- 
trix is given by the k x k symmetric matrix: 



iogC(x\e) 

ddidei 



(19) 



the average (...)= / £(X|0) . . . (i"X is taken over all 
possible data realizations, given the model parameters. 
One useful property of the Fisher information matrix is 
the Cramer-Rao theorem. It states that the variance 
for the parameter 9i about any unbiased estimator value 

cannot be inferior to (F~^)|/^, if the other parameters 
are estimated from the data as well and marginalized, or 
(Fji)"^/^, if all the other parameters are known. More- 
over, the information is additive, so that the total infor- 
mation allowed by two independent experiments is the 
sum of the information from each experiment separately. 

It is worth noting that the likelihood can be locally 
approximated by a multivariate Gaussian distribution. 
Taylor expanding to second order (the first non- vanishing 
term, up to an irrelevant additive constant) around its 

maximum 6 which, in the limit of large data sets, be- 
comes the best unbiased estimator of the actual param- 
eter set. Then the likelihood becomes: 



(20) 



exp 



-i(Xobs - X(0))^C(0)-i(Xobs - X(0)) 



where Xobs a-nd X(0) are respectively the observed and 
the theoretical values of the observables X and 









e=e- 



(21) 



is a positive semi-definite non-singular matrix, dubbed 
covariance matrix of the 9i. 

Therefore, under the hypothesis of a gaussian likeli- 
hood, the Fisher information matrix components are the 
expectation values of C^^{6). Accordingly, the inverse 
of the Fisher matrix is an estimate of the covariance ma- 
trix of the parameters C{6) « F^^. This shows that the 
Fisher matrix analysis represents only an approximation 
to the true likelihood and its results could significantly 
vary from the MCMC ones, especially if the parameters 
do have not a gaussian distr ibution. 

It can be shown (Tcgmar k et al.lll997| ) that, when the 
probability distribution is gaussian, the fisher informa- 
tion matrix (|19p reduces to 

F.yW-^Tr[C-iC,,C-iC„] , (22) 

if the means of the data are fixed. 

4. MCMC RESULTS 

4.1. Cosmological data 

In the case of the exponential potential (H^), likeli- 
hood contours and distributions for increasing (/3 > 0) 



or decreasing {[3 < 0) i/ mass are shown in Fig. [51 We re- 
port plots for Wj/, a and /3 which are the most significant 
parameters in our models, the constraints on the other 
parameters not considerably differing from the ACDM 
case. 

For both increasing and decreasing mass, no stringent 
constraints can be put on the coupling (3, whose upper 
limit exceeds |/3| = 5 yet at 68% confidence level (CL). 
On the contrary, the v density parameter is well con- 
strained at 95% CL: < 0.0034 (decreasing mass) and 
iOi, < 0.0032 (increasing mass). This limits are ^ 3 time s 
narrower than those found in iBrookfield et all (|2006b( ). 
In term of today neutrino mass these values attest an up- 
per limit to nil, — 0.32eV. The limits on a are centered 
around zero in both cases: —1.12 < a < 0.77 (decreasing 
mass) and —0.89 < a < 1.05 (increasing mass) at 1-a. 

Let us underline that the results show a common be- 
haviour in the /3 vs. plots, almost symmetric with 
respect to the sign of /3. In fact a degeneracy appears 
between the two parameters so that high are in agree- 
ment with low neutrino content, or equivalently low neu- 
trino mass. 

In Fig. [3] we reported likelihood contours for the case of 
the inverse power-law potential as in (I13p. This potential 



type w as not previously considered in IBrookfield et al] 
(|2006bf) . In the /? vs. lo^, plots the zero coupling is appar- 
ently excluded with statistical significance higher than 
95% confidence level (CL), with /? > 3 for the increas- 
ing neutrino mass and /? < —1.5 for the decreasing mass 
case at 2-cr. In particular, the growing neutrino mass case 
exhibits an explicit preferen c e for hi gh j3 values, in agree- 
ment with lAmendola et al.l (|2008af ) . For uji, only upper 
hmits can be found at 95% CL: < 0.0017 (decreasing 
mass) and a;^ < 0.0423 (increasing mass). 

4.2. External priors 

Complementary information to cosmological data can 
be obtained by taking into account priors on parame- 
ters coming from external measurements. In such a way 
a better knowledge of the other parameters can be ob- 
tained. 

In our analysis we fix the value of the neu- 
trino mass today (and as a consequence uj,^), us- 
ing two possible options, = 0.2 eV and TOj/ = 
0.3 eV . These values are within the range of 
the claimed i/g m ass detection in the H eidelberg- 
Moscow experiment (iKlapdor-Kleingrothaus et al. 200_J; 
iKlaEdor-Klcingrothaus 2005). Also the KATRIN experi- 
ment (jSturm..2011. ) is expected to constraint the value of 
the neutrino mass with a sensitivity in the sub-eV range. 

Beside fixing , we chose to fix the baryon density pa- 
rameter flbh^ ~ 0.022. as constraints on this parameter 
are not significantly modified by neutrino-DE coupling. 
Also the value for Hq = 72 km s~^Mpc~^ is kept con- 
stant. 

In the exponential potential case we are therefore left 
with a cosmological model requiring six parameters: /?, 
a, QcDAih'^, Zre, and ris. 

In Fig. m and [5] we show likelihood contours and dis- 
tributions for increasing and decreasing neutrino mass, 
respectively. The most significant result in both cases is 
the improvement in parameter constraints. In fact the 
limits on the coupling become more stringent, with an 
upper value |/3| ~ 2 3, and the distributions for a be- 



Fig. 2. — In the upper (lower) panel, likelihood contours and dis- 
tributions are shown for decreasing (increasing) neutrino mass in 
the exponential potential case. For the plots on the diagonal, dot- 
ted lines are mean likelihoods of samples, solid lines are marginal- 
ized probabilities. Similarly, for the coloured plots, black lines are 
1- and 2-cr contours of the marginalized probability distribution, 
while the colours refer to mean likelihood degradation from the 
top (dark red) to lower values. 



come narrower, but always compatible with zero at 68% 
CL. 

It is worth to notice that the colored plot in the upper 
panel in Fig. [4] shows a complex structure for the likeli- 
hood contours in which the upper value of the coupling 
still exceeds /3 = 5. This effect can be mainly ascribed 
to the peculiar features of the likelihood in Fig. [5] which 
are completely lost when considering a prior with higher 
mass. 

5. FISHER MATRIX FORECASTS 



Fig. 3. — As in Fig. [2] when the power-law potential is considerd. 

In the next few years new and more accurate datasets 
will be available. In this section we show Fisher ma- 
trix analysis results from the combination of the CMB 
anisotropics and the tomographic weak lensing (TWL) 
spectra. 

We considered the case of the exponential potential 
with an exponential coupling. The fiducial parameters 
6a are shown in Tab. [TJ This two sets correspond to 
the maximum area of the likelihood as determined in 
par. 14.11 and are chosen to be exactly the same, with the 
only exception of the sign of /3. 

T he Fisher matrix for CMB power spectrum is given 
bv (iZaldarriaga fc Selial3 119971 : iZaldarriaga et all 119971 : 
iRassat et al.ll2008[ ): 



pCMB. = ^ ^ dCx,i 



I X,Y 



.cov^V^ 



(23) 



6 






Fig. 4. — Upper (lower) panel shows posterior constraints for Fig. 5. — Upper (lower) panel shows posterior constraints for 

the exponential potential with increasing v mass. Here we fixed the exponential potential with decreasing v mass. Here we fixed 

mv = 0.2eV (m^ = 0.3eV), (.j^ = 0.022, Ho = 72 = 0.2eV (m^ = 0.3eV), (.j;, = 0.022, Ho = 72 



TABLE 1 

Fiducial cosmological parameters for the exponential 
potential, consistent with fig. [2] 



Parameter 




^CDM 




ns 


In As 


Value 


0.0224 


0.1097 


69 11 


0.95 


2.09x10-^ 


Parameter 




a 








Value 


0.0010 


0.02 


±1.9 







where Cx,i is the harmonic power spectrum for the 
temperature-temperature (X = TT), temperature- 
E-polarization {X = TE) and the E-polarization-E- 
polarization [X = EE) power spectrum. Under the as- 
sumption of a Gaussian hkehhood distribution, the co- 



variance COVjj-y of the errors of the various power spec- 
tra is given in terms of Cx,i with 



COVt,t = /; (Ct,( + WT^B-'^f 
COY E.E = fi {C'Ej + Wp'B-^y 



COYte,te ~ fi 



TE,, 



(24) 
(25) 
(26) 



I I V^B,i 
COVt.B = flCTE,l 

CO\t,te = fiCTE^i {Ct,i + W^^BY"") 
COY E,TE^flCTE,l {Ce,i + Wp'B-^) 



(27) 
(28) 
(29) 
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where /; = (2i+i)/sky ^^'^ ^t.p = (o'T.p^'fwhm) ^ is the 
weight per soUd angle for temperature and polarization, 
with a l-cr sensitivity per pixel of (Tt,p with a beam 
of ^twhm extend. The beam window function is given 
in terms of the full width half maximum (fwhm) beam 
width by Bi = exp + I)ef„^^j2ln256) and /,ky is 
the sky fraction. 

We will use the for t hcomi ng Planck mission 
(|The Planck CollaborationI I2006f ) specifications for 
the measurements of the CMB temperature and po- 
larization spectra. Note that Eq. usually includes 
a summation over the Planck frequency channels. 
However we conservatively assume that we will only use 
the 143 GHz channel as science channel, with the other 
frequencies used for foreground removal (not treated in 
this paper). This channel has a beam of 0fwhm = 7.1' and 
sensitivities of ar = 2.2^K/K and ap = A.2ijlK/K. To 
account for the galactic plane cut, we take /sky = 0.80. 
Note we use as a minimum £-mode, ^min = 30 in order 
to avoid problems with polarization foregrounds and 
subtleties for the modeling of the integrated Sachs- Wolfe 
effect. 

The Fisher matrix for TWL reads: 

fZ"- = (30) 

{21 + \)M dP,j,i 1 dPk^,e 1 
hky 2 96»„ ^■'^■^ dBfi ■ 

where Pij^i are the components of the non-linear weak 
lensing power s pectrum for the i-th and j-th bin 
(|Hu fc Jainll200l . In the expression (|30p a summation 
over repeated indices is implicit. The redshift bins have 
been chosen such that each contains the same amount of 
galaxies. We consider the 5 bin case as reference case. 

The n on-linear correctio ns are calculated using 
HALOFIT (jSmith et al.l l2003f ). This procedure is suit- 
ably fitted to ACDM N-body simulations. Therefore 
it could lead to errors of the order of 20% on Fisher 
outp uts, if used for models different from ACDM (see 
e.g. iCasarini et al.l ()2011[ )). However, in the absence 
of suitable extensions for MaVaN, we assume HALOFIT 
as the procedure for non-linear corrections, reporting re- 
sults up to Iraax = 1000, in the mildly non-linear regime. 
This conservative choice also prevent us from considering 
a regime in which baryons could strongly affect matter 
pow er spectra and, as a consequence, weak lensing spec- 
tra (|Casarini et al.ll2012D . 

The TWL survey parameters are fixed accord- 
ing to the forthcorn ing Euclid missioifl specifications 
(jLaureiis et al.l [201 It ). The survey area covered by the 
experiment is 15000 deg^, while the density is 30 galax- 
ies per arcmin^. The distribution of the galaxy number 
on the redshift and solid angle is n(z) = n^z^ e^-^ I 
with a median redshift z = 0.9. The photometric redshift 
error is 0.05 (1 4- z). 

The error estimation of DE and neutrino density pa- 
rameters is reported for a MaVaN model with an expo- 
nential coupling for an increasing and decreasing neu- 
trino mass. As shown in Tab. [2l TWL is able to put 
constraints stronger than CMB on /? and oj^. with an 

^ http://www.euclid-ec.org 



TABLE 2 

l-cr ERROR ESTIMATIONS OF COSMOLOGICAL PARAMETERS FOR AN 
EXPONENTIAL POTENTIAL WITH /3 < 0, USING PlANCK-LIKE 

CMB AND Euclid-like Tomographic Weak Lensing (TWL) 

DATA WITH Iraax = 1000, ALONE OR IN COMBINATION. 











TWL 


0.00043 


0.024 


1.17 


CMB 


0.00064 


0.014 


1.28 


TWL+CMB 


0.00025 


0.012 


0.69 




TABLE 3 






The same as 


Tab [2] but 


WITH 


/3 > 0. 


Dataset 








TWL 


0.00045 


0.028 


2.04 


CMB 


0.00054 


0.019 


1.71 


TWL+CMB 


0.00026 


0.012 


1.15 



improvement of -^10% and >30%; respectively. On the 
contrary CMB is more efficient in constraining a, gaining 
^50% over TWL. The combination of the two observ- 
ables can improve constraints on all parameters. In this 
case, the main result is for the estimate of constraint 
to 0.00025, (Tm„ ^ 0.02eV. Analogous comments can be 
made for the /3 > in Tab. [31 with the only difference 
that CMB can constrain /3 better than TWL. 

In Fig. [6l we show l-cr and 2-cr likelihood contours for 
MaVan's model with an exponential coupling for both an 
increasing and decreasing neutrino mass. Here we report 
results for Wi/, /3 and a, after marginalizing over the other 
parameters. It is quite clear that neither Planck nor 
Euclid alone will be able to constrain a non-zero coupling 
/3 ^ C(l), only TWL being slightly more efficient than 
CMB in the /3 < case. In this case the combination of 
the two can exclude a null coupling with a significance 
higher than 2-cr. 

Let us note that we performed a conservative analysis 
of TWL Fisher matrix, using multipoles up to Imax = 
1000. The capability of TWL in constraining /3 and the 
other parameters will be significantly increased by the 
use of higher multipoles, delving into the deep non-linear 
regime. For instance Fig. [7| presents the confidence el- 
lipses obtained including multipoles up to imax = 5000 
in the TWL Fisher matrix. Estimated errors are clearly 
reduced of a factor ~2 with respect to Fig. [7| for Weak 
Lensing and, as a consequence, for the combined case. 
This figure shows that Euclid would be able to exclude 
the zero value both for /3 and uju with high statistical sig- 
nificance. However in this case an appropriate extension 
of HALOFIT would be necessary, including both appropri- 
ate non-linear correction for MaVaN models and baryon 
physics. 

Despite this limitation, a very interesting results is the 
fact that TWL, alone or in combination with CMB, is 
able to find a lower value to VL^lr? with an error half of 
the CMB error if taken alone. Moreover a strong degen- 
eracy can be found in the a;,y-/3 plane for what concerns 
the TWL ellipses, also if considered in combination with 
CMB. This effect appears in both cases in a symmetric 
way, showing that lower |/3| values allow higher neutrino 
mass. The Vly\T?-i5 degeneracy could be eventually bro- 



-0.02 0.020.040.06 



-0.02 0.020.040.06 




^-0.02 0.020.040.06 ^0 1 „ 2 




Fig. 6. — Upper (lower) panel shows likelihood contour forecast 
at 68% and 95% confidence levels for the exponential potential, 
in the decreasing (increasing) mass case. We marginalized over 
the other parameters. The green lines are obtained considering 
a tomographic weak lensing Euclid-like survey, with 5 bins and 
(■max = 1000. The blue lines are for a Planck-like CMB mission. 
The red lines are for a combination of the two. The crosses indicate 
the fiducial values 

ken from external priors on rriu- On the other hand, no 
noteworthy remark can be said about the a parameter, 
whose constrains remains compatible with zero also for 
the combination of the two observables. 

6. DISCUSSION AND CONCLUSIONS 

In this work we focused on the hypothesis that the 
origin of the cosmic acceleration can be attributed to a 
quintessence scalar field coupled to massive neutrinos. 
Firstly, we updated and extended parameter constraints 
using the most recent available data from SNIa, CMB 
and LSS observations. We considered both an exponen- 




FiG. 7. — As FigElwith 5 bins and i^ax = 5000. 

tial and a power law potential, with increasing (/3 > 0) 
or decreasing (/? < 0) neutrino mass. 

In none of these cases the cosmological data placed 
strong constraints on the coupling parameter in the low 
coupling range. Neither for the neutrino density parame- 
ter, on which only upper limits can be placed. The main 
outcome of the analysis was that /3 values ^ 0{1) are 
compatible with actual data, together with a neutrino 
mass rrii, < 0.32eV. 

Therefore at the moment we do have not enough in- 
formation to exclude a possible coupling between neu- 
trino and DE field, neither in the low or in the high 
coupling regime. New and precise data from observables 
related to the recent evolution of the universe are neces- 
sary together with a deeper understanding of the MaVaN 
theory dynamics. In this sense the forthcoming Euclid 
weak lensing and galaxy clustering data will represent 
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the turning point for the future progress in cosmology. 

With this aim, a Fisher matrix study was accom- 
phshed, considering future data release from missions 
such as Planck for CMB spectra and Euclid for tomo- 
graphic weak lensing spectra. The latter experiment 
will be more efficient in improving present constraints 
on LOy, even considering a cautious use of multipoles up 
to tmax = 1000. This choice prevent us from including 
highly non-linear features, which are not correctly pre- 
dicted by HALOFIT for MaVaN theories. Combining Eu- 
clid data with complementary information from Planck 
will exclude a zero coupling at more than 2-cr. It is worth 
to mention the strong degeneracy between /3 and ujy, 
which could eventually be broken by an external prior 
on the neutrino mass. 

A crucial point for the forthcoming Euclid mission will 
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